An algorithm for computing the integral closure 
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In this paper we describe an algorithm for computing the normalization for certain rings. 
This quite general algorithm is essentially due to Grauert and Remmert 0, [§, and seems 
to have escaped the attention of the computer algebra specialists until now. As I am not a 
computer algebra specialist myself, I do not know whether this algorithm is fast. Grauert 
t^. ! and Remmert proved a normality criterium in order to give a simple and nice proof of a 
theorem of Oka, who proved that the non-normal points of an analytic space is an analytic 
space itself. We just reformulate their results, so that it better suits our purposes. 

u 
< 



To fix the notation, let R be a Noetherian and reduced ring, R the integral closure (also 
called normalization) of R. We consider the set: 

NNL := {p G Spec(.R) : R p is not normal} 



Here NNL stands for non-normal locus. Let / be an ideal of R containing a nonzerodivisor. 
■ We have canonical inclusions: 

O 

g ; R C Hom R (/, I) C R 

t> . 

Q\ . The first inclusion is the map which sends an element of R to multiplication with this element. 

The second inclusion is sending G Honis(I, /) to ^p- for any element f £ I which is a 

nonzerodivisor of R. It is easily checked that the map is independent of the choice of /. 

That we in fact land in R can be found in any textbook which has included integral closure 

i as a topic. 

JSC; 

Theorem 1. /TJ/ pp. 220-221, /fj/, pp. 125-127. Assume that the ideal I contains a nonze- 
£> . rodivisor, and has the following property: 

^ ■ NNL C V(I) 

where V(I) = {p E Spec(-R) : / C p} denotes, as usual, the zero set of I . Suppose moreover 
that I has the property 

Rom R (I, I) = Hom iJ (J, R) n R (*) 

Then one has the following normality criterium: 

R = Homft(/, J) <^=^> R is normal 

Proof. The implication is trivial. To prove the converse, let h — jj G R. Consider the 
following ideal in R 

{geR-.hgeR} 



1 



Its zero set is called the "pole set" of h: 



P(h) := {p G Spec(R) :h<£R p } 

It is immediate that P(h) C NNL. Let J be the ideal of P(h). There exists a c > such 
that hJ c C -R, by the Nullstellensatz. By the Nullstellensatz again y/l C J. Therefore there 
exists a d > such that hl d C -R. Let d be minimal with this property. We claim d — 1. 
Suppose the converse, i.e. d > 1. Then there exists an a G with ft,a ^ -R. Furthermore 
ha E R and (ha) I C -R. By assumption: 

i2 = {/i G -R : W C J} = {h G -R : hi C /} 

so that (ha) G J. Therefore ha & R after all, a contradiction. □ 

We have to find an ideal which satisfies condition (*). This is provided by: 

Theorem 2. Every radical ideal I containing a nonzerodivisor safisfies condition (*). 

Proof. The proof is in and ||, but because it is so nice and simple we give it here. Let 
h G -R, so we have an equation: 

h n = ao + a\h + . . . + a n _i/i n ~ ; a« G -R 
If hi C -R, then we have for all / G /: 

(W = a r + aih + ■ ■ ■ + a n ^h n - l f G / 
As / is supposed to be reduced it follows that hf G /, and that is what we had to prove. □ 
These results give rise to the following algorithm: 

ALGORITHM 

INPUT: A reduced noetherian ring R. 

OUTPUT: The normalization R of R. 

STEP 1: Determine an ideal I with NNL CV(I). 

STEP 2: Compute the radical \fl of I. Put I :— \fl. 

STEP 3: Take an / G I, and compute J := Ann(f). If J = 0, GOTO STEP 5. 
STEP 4: Put R := R/(f) © R/J and GOTO STEP 1. 

STEP 5: Compute Hom fl (J, I). If R = Hom fl (J, I) then put R := R and STOP. 
STEP 6: Set R := Rom R (I,I) and GOTO STEP 1. 

This algorithm stops exactly when the normalization R is finitely generated as an R- 
module, so for example for affine rings, due to a classical result of E. Noether. 

Some remarks are in order. 

1. To determine an ideal / with NNL C V(I) one can take any I which contains the 
non-regular locus of -R. 
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2. Algorithms for computing the radical of an ideal are described in and probably will 
take the longest in this algorithm. 

3. In step 4, R C R/(f) © Rj J is indeed an inclusion, as (/) fl J = follows from the 
assumption R is reduced. This extension is also finite. To see this, take a g 6 J with 
f + g a nonzerodivisor (prime avoidance). The extension: 

R C © = R[X]/(X 2 - X, X(f + g)-f) 

obviously is finite. As R C R/(f) © Rj J is a quotient of R/(f) © R/(g) (in fact it is 
equal) it is a finite extension too. 

4. I think it should be possible to extend this algorithm to an algorithm which computes 
the primary decomposition for a radical ideal. 

We finish with describing the ring structure of Homjj(7,7), which essentially is due to 
Catanese, see . Take generators uq := 1, m, . . . , u t of Hohir(J, I) as i?-module. Consider 
the map: 

R ■ X © R ■ X 1 © . . . © R ■ X t Hom R (/, J) 



Xi H-> Wj 

Computing the kernel of the map gives "linear equations": 

t 

Li = a ijXj = a^j G R; i — 1, . . . , s 

3=0 

Because Hom^(J, I) is a ring, we have that u^j for all 1 < i < j < t is in Hom^(J, J) again. 
We therefore can find elements ftijk G -R such that: 

k=0 

giving us tt£±ll "quadratic equations": 

t 

Qij := XiXj — PijkXk 

k=0 

For the easy proof of the following theorem we refer to Catanese [|l| . 

Theorem 3. Put X = 1, and consider the ideal J C R[Xi, . . . , X t ] generated by the Li, i = 
1, . . . , s and the for 1 < i < j <t. Then we have a ring isomorphism: 

Rom R (I,I) = R[X u ... ,X t ]/J 



Acknowledgment I am very grateful to Wolfram Decker for bringing to my attention 
that this algorithm was unknown to the computer algebra specialists. 
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